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В рамках приближения локальной плотности выполнены аЬ іпіііо расчеты электронной структуры чистого Рб и системы РбН х 
(х=1, 2, 3). Вычислена полная энергия системы РсІН х , для случаев различной координации атомов водорода (окта- и тетраэдри- 
ческой), сделан вывод о наиболее вероятном их местоположении в решетке металла. В приближении постоянного матричного 
элемента вычислена мнимая часть функции диэлектрической проницаемости е 2 (ю). Обнаружено, что растворение водорода в 
палладии увеличивает значения функции е 2 (со) в исследованной области энергий от 2 до 24 эВ. Поэтому в случае радиационно- 
го воздействия на систему РбН х можно ожидать интенсивного возбуждения электронной подсистемы кристалла и, как следствие, 
снижения потенциальных барьеров для перемещения атомов водорода. 


Введение 

Уникальная способность палладия растворять 
большое количество водорода была открыта Т. Грэ- 
хемом еще в 1886 г. Растворение водорода значи- 
тельно изменяет физические свойства палладия 
[1,2]: увеличивает параметр решетки, снижает элек- 
тропроводность и магнитную восприимчивость, по- 
вышает твердость и прочность, снижает пластич- 
ность, приводит к охрупчиванию. Система РсІН, яв- 
ляется диамагнетиком и обладает сверхпроводимо- 
стью, тогда как чистый Рё - парамагнетик и, как из- 
вестно, не проявляет сверхпроводящих свойств [3]. 
Очевидно, что эти изменения свойств палладия свя- 
заны с модификацией его атомной и электронной 
структуры при растворении в нем водорода. 

Растворяя водород, палладий сохраняет тип своей 
кристаллической решётки. При малых концентра- 
циях водорода (х<0,03) образуется так называемая сс- 
фаза РсІН,. При высоких концентрациях водорода 
(х>0,6) раствор переходит в /1-фазу с характерным для 
фазового перехода первого рода скачкообразным из- 
менением постоянной решетки. Растворяясь в палла- 
дии, атомы водорода располагаются в междоузлиях 
его гранецентрированной кубической (ГЦК) кри- 
сталлической решетки. Как известно, ГЦК решетка 
имеет два типа междоузлий: тетраэдрические и окта- 
эдрические. Основываясь на данных по дифракции 
нейтронов [4-6], был сделан вывод о том, что атомы 
водорода занимают октаэдрические междоузлия в 
ГЦК решетке палладия, образуя моногидрид РсІН со 
структурой ІЧаСІ. Однако имеются данные о концен- 
трационной и температурной зависимости электри- 
ческого сопротивления моногидрида РёН, а также 
нейтронно-структурные исследования [3, 6], свиде- 
тельствующие о переходе атомов водорода из окта- в 
тетрапоры при температурах порядка 50 К. 

В ходе недавних экспериментов получена инфор- 
мация о поведении водорода в решетке палладия под 
действием ионизирующего излучения. Так, получены 
экспериментальные данные о радиационно-стиму- 
лированной миграции и десорбции водорода из па- 
лладия, в том числе и в атомарном (неравновесном) 
виде в широком диапазоне температур, под действи- 
ем электронного, рентгеновского и у-облучения 


[7, 8]. По мнению авторов этих работ, такое поведе- 
ние водорода обусловлено его способностью погло- 
щать и аккумулировать энергию ионизирующего из- 
лучения на время, достаточное для совершения нес- 
кольких скачков из своего первоначального положе- 
ния. Однако никаких микроскопических механизмов 
этого явления до настоящего времени не предложе- 
но. Ставя перед собой задачу создания такой микро- 
скопической модели, на первом этапе мы попыта- 
лись выяснить, как влияет присутствие водорода в 
решетке металла на способность электронной подси- 
стемы кристалла поглощать энергию внешнего ра- 
диационного воздействия. Для ответа на поставлен- 
ный вопрос необходима информация об электрон- 
ной и атомной структуре исследуемого гидрида. Зная 
закон дисперсии Е(к) для элементарных возбужде- 
ний твёрдого тела (в данном случае это система элек- 
тронов), можно вычислить функцию отклика систе- 
мы на внешнее возмущение (частотную зависимость 
комплексной диэлектрической проницаемости). 

В настоящей работе на основе расчетов аЪ іпіііо 
полной энергии исследовано образование соедине- 
ний РсІН, для х=1, 2, 3 из элементных компонен- 
тов: Рё и Н 2 . Рассчитаны плотности электронных 
состояний ЩЕ), дисперсионные кривые Е(к), ра- 
спределения зарядовой плотности и мнимая часть 
диэлектрической функции е 2 (со). 

Метод и детали расчета 

В рамках теории функционала плотности ли- 
неаризованным методом присоединенных плоских 
волн [9, 10] проведены самосогласованные расчеты 
электронной структуры системы РёН х и чистого 
Рё. Обменно-корреляционные эффекты рассма- 
тривались с использованием обобщенного гради- 
ентного приближения РЛѴ91 [11]. Радиусы 
МТ- сфер атомов Рё и Н были выбраны равными 
2,3 и 1,0 ат. ед. соответственно, что позволяло ато- 
мам водорода свободно помещаться в тетра- и ок- 
таэдрическом междоузлиях ГЦК решетки палла- 
дия. В разложении волновой функции внутри 
МТ-сфер по сферическим гармоникам учитыва- 
лись члены вплоть до /=8 для Рё и до 1=4 для Н. 
При разложении волновой функции по присоеди- 
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ненным плоским волнам учитывалось до 120 ба- 
зисных функций на атом. На каждой итерации са- 
мосогласованна собственные значения гамильто- 
ниана рассчитывались в 40 ^-точках неприводимой 
части зоны Бриллюэна (НЧЗБ). Для построения 
плотности электронных состояний ЩЕ) проводи- 
лось интегрирование в НЧЗБ линейным методом 
тетраэдронов на сетке из 90 Л-точек. Сглаживание 
функции ЩЕ) осуществлялось гауссианами со 
среднеквадратичной шириной 0,14 эВ. Расчет мни- 
мой части диэлектрической функции е г (со) прово- 
дился в приближении постоянного матричного эл- 
емента на сетке из 1200 Л-точек в НЧЗБ. 

Результаты и обсуждение 

Процесс растворения водорода в палладии 
можно представить следующим образом: 

тРё + Н 2 <-> тРёН х - А Е, 

где теплота растворения водорода 

А Е = ^(РёН х )- Я(Рф- |я(Н 2 ). 

Здесь ДРёНД и ЯРё) - полные энергии систе- 
мы РёН х и чистого Рсі, соответственно. Я Н) - 
энергия, необходимая для диссоциации молекулы 
Н 2 . Полученное значение для энергии диссоци- 
ации ДН 2 ) составило -30,96144 эВ. В табл. 1 при- 
ведены результаты расчета параметров решетки и 
теплоты растворения рассматриваемых гидридов. 
Видно, что результаты, полученные для параметров 
решетки, хорошо согласуются как с эксперимен- 
тальными [2, 14, 15], так и с расчетными данными 
других авторов [12, 13, 15]. Что касается теплоты 
образования гидридов, то здесь ситуация оказалась 
непростой. Так, если не учитывать энергию нуле- 
вых колебаний атомов водорода в металле, энерге- 
тически наиболее выгодной будет тетраэдрическая 
их координация. Учет энергии этих колебаний в 
гармоническом приближении приводит к тому, что 
энергетически более выгодным для атомов водоро- 
да становится октаэдрическое междоузлие. 


Таблица 1. Параметры решётки и теплота растворения ги- 
дридов палладия 


Си- 

сте- 

ма 

Параметр решётки, А 

Теплота растворения водоро- 
да, эВ 

На- 

стоя- 

щая 

работа 

Теория 

Экспе- 

римент 

Настоя- 
щая ра- 
бота 

Теория 

Экспери- 

мент 

Рсі 

3,93 

3,97 [12] 

3,89 [2] 

- 

- 

- 

РсІН 

4,08 

4,06 [13] 
4,15 [12] 

4,09 [2] 

-0,533 

-0,462* 

-0,143 [13] 
-0,500 [15] 

-0,21 [14] 
-0,122 [15] 

РсН 

4,47 

4,40 [13] 

- 

-0,692 

-0,238* 

-0,216 [13] 

- 

РсІНз 

4,66 

4,64 [13] 

- 

1,311 

- 

- 


Звёздочкой отмечены результаты, полученные с учётом энер- 
гии нулевых колебаний атомов водорода в гармоническом 
приближении 


Таблица 2. Локальное орбитальное распределение валент- 
ного заряда 0 (в ед. заряда электрона) для чисто- 
го палладия и систем РбН х при х=1, 2, 3. Значения 
плотности электронных состояний на уровне Фер- 
ми Л/(2У (в ед, сост./эВ/эл. яч.) 


Си- 

стема 

Атом 

0, 

0 Р 

Оз 

Омт 

Оіш 

Ош 

М Е г ) 

рб 

Рсі 

0,274 

0,201 

7,796 

8,299 

1,701 

10 

2,64 

РсІН 

Рсі 

0,214 

0,210 

7,908 

8,356 

2,261 

11 

0,49 

Н (окт.) 

0,363 

0,018 

0,001 

0,383 

1 1 

РбН 2 

Рсі 

Н (тет.) 

0,197 

0,392 

0,236 

0,017 

7,907 

0,001 

8,364 

0,411 

2,814 

12 

0,13 


Рсі 

0,210 

0,198 

7,892 

8,318 




РсІНз 

Н (окт.) 

0,305 

0,010 

0 

0,315 

3,615 

13 

1,42 


Н (тет.) 

0,361 

0,014 

0,001 

0,376 





В табл. 2 использованы следующие обозначе- 
ния: 0 Ш - полный заряд в элементарной ячейке, 
Ош - заряд в межатомной области, 0 МТ . О,. 0 Р (2, - 
заряд в МТ-сферах и его 5-, р- и «'-компоненты. 
ЩЕ Р ) - значение плотности электронных состоя- 
ний на уровне Ферми. 

На рис. 1 и в табл. 2 представлены рассчитан- 
ные в работе плотности электронных состояний 
для чистого Рё и систем РёН, (х=1, 2, 3) с различ- 
ной координацией атомов водорода. Полученные 
результаты хорошо согласуются с данными, полу- 
ченными ранее другими авторами [16, 17]. 

Из сравнения зонных структур чистого Рё и си- 
стем Рё-Н можно видеть характерные особенно- 
сти электронной структуры металлического гидри- 
да, возникающие в результате растворения водоро- 
да: 

1. заполнение «1-зоны Рё за счёт электронов ра- 
створённого водорода; 

2. понижение дна зоны проводимости за счет ги- 
бридизации 5-состояний водорода с 5-, р- и 
«/-состояниями атомов Рё; 

3. сдвиг уровня Ферми вверх относительно потол- 
ка «'-зоны Рё, что приводит к уменьшению 
плотности электронных состояний системы на 
уровне Ферми. 

Из рис. 1 также видно, что появление водорода 
в решетке Рё приводит к появлению псевдощели в 
электронном спектре системы Рё-Н. Для соедине- 
ния РёН эта псевдощель лежит в области энергий 
~4 эВ. С ростом концентрации водорода она опу- 
скается вниз по энергии за счет заполнения элек- 
тронных состояний электронами атомов водорода. 
Для системы РёН 2 с тетраэдрической координаци- 
ей атомов водорода уровень Ферми попадает в 
псевдощель, и соединение приобретает полуметал- 
лические свойства. Дальнейший рост концентра- 
ции водорода приводит к тому, что система РёН 3 
снова становится металлом, т.к. псевдощель опу- 
скается ниже уровня Ферми более чем на 1 эВ. На- 
личие псевдощели в электронном спектре системы 
Рё-Н говорит о появлении ковалентной соста- 
вляющей сил связи между атомами палладия и во- 
дорода в дополнение к металлическому типу связи, 
характерному для чистого палладия. 
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Е - Ер, )В Е - Е р , эВ 

Рис. 1 . Плотность электронных состояний чистого РА и соединений РАН , 

Известно, что уменьшение плотности состоя- чистом палладии (табл. 2), поэтому для них следует 
ний на уровне Ферми приводит к снижению экра- ожидать увеличения радиуса экранирования и, как 
нирующей способности электронов зоны проводи- следствие, увеличение времени жизни электрон- 
мости. Действительно, согласно модели Фриделя ных возбуждений. 

[2] избыточный заряд в металле экранируется элек- На р ИС 2 приведены дисперсионные кривые 
тронами зоны проводимости, плотность которых дд вдоль симметричных направлений неприво- 
экспоненциально спадает с увеличением расстоя- димой части зоны Бриллюэна, рассчитанные в 
ния и характеризуется радиусом экранирования рамках теории функционала плотности линеаризо- 
(длиной спада) ~[./Ѵ(іѵ)]“ 1/2 . Так, если плотность со- ванным методом присоединенных плоских волн, 
стояний велика — радиус экранирования как ив Из р ИС у Нка видно, что растворение водорода в па- 
чистом палладии мал, поскольку даже ближайшие дладии не приводит к появлению дополнительных 
соседи не чувствуют замены атома матрицы атомом энергетических зон, но заметно изменяет их ход и 
водорода [2]. Для рассматриваемых систем Р<Ш Х взаимное расположение. За счет заполнения ды- 
плотность состояний на уровне Ферми ниже, чем в рочных состояний в А - зоне палладия электронами. 
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Рис. 2. Дисперсионные кривые для Рб и соединений РбН г 


отданными водородом в общую зону проводимости 
системы Рсі — Н, г/-зона становится более узкой и 
опускается по энергии. 



Энергия, эВ 

Рис. 3. Мнимая часть диэлектрической функции Еі(со) для 
чистого Рб и систем РбН х 


На рис. 3 приведена рассчитанная мнимая часть 
диэлектрической функции 8 2 (со). Эта величина 
устанавливает связь между структурой зон и опти- 
ческими константами вещества. В частности, она 
характеризует способность вещества поглощать 
электромагнитное излучение. 

В литературе не удалось обнаружить каких-ли- 
бо данных о рассчитанной нами величие е 2 {со) для 


системы Рё-Н. Сравнение же результатов для чи- 
стого палладия с данными полуэмпирических рас- 
чётов [17, 19] и экспериментальных работ [20] по- 
казало неплохое согласие. Полученная зависи- 
мость е 2 (а>) в целом отражает ход эксперименталь- 
ной кривой [20], но имеет отличие от результата 
расчётов [17, 19]. Так, в полученных нами результа- 
тах отсутствует пик в области энергий порядка 
4,5 эВ. Отметим, что в эксперименте этот пик так- 
же не наблюдается. Из рис. 3 видно, что значения 
диэлектрической функции е 2 (со) рассматриваемых 
соединений РёН х превышают соответствующие 
значения е 2 (со) для чистого Рё в широком диапазо- 
не энергий (примерно до 24 эВ). Таким образом, 
можно говорить, что при растворении водорода в 
палладии повышается способность системы Рсі-Н 
поглощать энергию электромагнитного излучения. 
Это обстоятельство совместно с уже отмеченным 
выше увеличением времени жизни электронных 
возбуждений в системах Рё-Н может существенно 
изменить характер поведения соединений РёН х в 
поле ионизирующего излучения. 

Известно, что если в кристалле или его фраг- 
менте энергия электронов повышена, можно ожи- 
дать снижения потенциальных барьеров для мигра- 
ции атомов [21]. Так, эффективно поглощая энер- 
гию радиационного воздействия, электронная под- 
система соединений РёН* перейдёт в возбуждённое 
состояние. Это изменит потенциальных барьеры 
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для миграции атомов и, в первую очередь, для наи- 
более легких. Если время жизни этих возбуждений 
будет достаточно большим, атомы успеют совер- 
шить несколько скачков, прежде, чем возбуждение 
релаксирует. Увеличением продолжительности ра- 
диационного воздействия можно добиться того, 
что атомы водорода начнут интенсивно переме- 
щаться по объёму кристалла и выходить за его пре- 
делы. 

Заключение 

Проведены первопринципные расчеты элек- 
тронной структуры чистого палладия и системы 
Рё-Н. Определены равновесные значения параме- 
тров решетки ряда соединений РёН х для х=1, 2, 3 
при различных координациях и концентрации во- 
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дорода в решетке палладия. Получено хорошее со- 
гласие рассчитанных величин с экспериментальны- 
ми и теоретическими данными других авторов. 
Проведен расчет мнимой части диэлектрической 
функции системы РёН х . Из анализа полученных ре- 
зультатов сделан вывод, что при растворении водо- 
рода в палладии увеличивается способность элек- 
тронной подсистемы металла поглощать энергию 
внешнего электромагнитного воздействия и сохра- 
нять её в течении более длительного времени, чем в 
чистом палладии. Выдвинуто предположение, что 
действие ионизирующего излучения на систему 
Рсі-Н может приводить к миграции атомов водоро- 
да и к их выходу из металла в неравновесном (ато- 
марном) состоянии за счёт снижения потенциаль- 
ных барьеров, обусловленного возмущением равно- 
весной зарядовой плотности в объёме образца. 
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